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Procede de synchronisation fine sur un signal r gu 
d'un canal d transmission 
La presente invention concerne un proc6d6 de 
synchronisation fine sur un signal de reception 
5 correspondant & un signal de reference £mis dans un canal de 
transmission . 

Dans un systeme de transmission, notamment par ondes 
radios, un emetteur emet un signal de reference dans un 
canal de transmission & destination d'un recepteur. Une d s 

10 premieres operations que doit r6aliser le recepteur est la 
synchronisation sur le signal de reception. Ce probleme est 
bien connu de l'homme du metier et il n'est done pas 
necessaire de rappeler ici les diff^rentes techniques qui 
sont employees pour obtenir cette synchronisation. 

15 Lorsque les signaux sont num6riques, on a coutum 

d'6valuer I'ecart de synchronisation au moyen d'une unite d 
temps, le temps bit, qui est l'6cart temporel s6parant deux 
bits successifs d'un signal. Or, il apparalt que les 
solutions disponibles dans I'etat de I'art ne permettent pas 

20 d'acqu^rir une synchronisation beaucoup plus precise que le 
temps bit. 

Cette precision peut s'av^rer insuffisante dans 
certains cas. En effet, e'est la synchronisation qui donne 
le temps de trajet du signal dans le canal de transmission 

25 ou, autrement dit, le temps de transmission entre 1' Emetteur 
et le r6cepteur. Cette donn6e est importante car, dans un 
systfeme de radiocommunication duplex od une station de bas 
est en communication avec un terminal, chaque 6quipement 
6tant pourvu d'un €metteur-recepteur , le terminal doit 

30 fonctionner de sorte que le signal qu'il 6met arrive a un 
instant precis en reference a la base de temps de la station 
de base. Pour ce faire, il faut naturellement que le 
terminal connaisse le temps que met ce signal pour parvenir 
a la station de base. 

35 D' autre part, c temps de transmission reflfet 

directement la distance separant l'6m tt ur du rec pt ur. On 
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comprend bien que la precision de cette distance est 
fondamentale lorsqu^l s'agit de proceder a la localisation 
du terminal en reperant sa position relative par rapport & 
une ou plusieurs stations de base. Le probleme general de la 
5 localisation d f un terminal est une preoccupation trfes 
actuelle du fait de ses applications au nombre desquelles on 
citera, par exemple, les strategies de changements de 
cellules ("handover") dans les reseaux cellulaires. On 
mentionnera egalement le domaine de la s6curit6, qu'il 
10 faille situer g6ographiquement la provenance d'un appel 
d'urgence ou bien la position d'un vehicule vol6 6quip6 du 
terminal . 

La pr6sente invention a ainsi pour objet un proc6d6 de 
synchronisation dont la precision est bien superieure au 
15 temps bit. 

Selon l 1 invention, le proc€d€ de synchronisation fine 
sur un signal de reception correspondant & un signal de 
reference emis dans un canal de transmission comprend les 
Stapes suivantes : 
20 - selection d'un signal source produisant un signal de 

caract6risation suite & son passage dans le canal de 
transmission , 

- etablissement d'une matrice de caracterisation pour 
estimer la covariance du signal de caracterisation, 

25 - identification des valeurs propres dominantes qui 

sont les valeurs propres les plus eievees de cette matric 
de caracterisation, 

- calcul de la fonction de correlation du signal 
source avec la somme des vecteurs propres associes aux 

30 valeurs propres dominantes, 

- recherche du premier maximum de cette fonction de 
correlation . 

Lorsque l 1 increment de temps adopte pour le calcul de 
la fonction d correlation est choisi suffisamment petit, en 
35 tout cas bien inf6rieur & un temps bit, c proc6d6 permet 
d'obtenir un tr^s bonne precision. 
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S Ion une premier option, le nombre de ces valeurs 
propres dominantes est predetermine, Typiquement , ce nombre 
represente de l'ordre de 20 a 30 % de la dimension de la 
irtatrice de caract6risation. 
5 Selon line deuxidme option, le rapport de la somme des 

valeurs propres dominantes h la somme de toutes les valeurs 
propres est sup6rieur ou egal a un nombre predetermine. Ici, 
le nombre retenu sera souvent superieur a 90 %, 95 % par 
exemple . 

10 Selon une troisifeme option, le proced<§ comprenant d 

plus une 6tape d' estimation du bruit additif dans le canal 
de transmission, les valeurs propres dominantes sont telles 
que leur somme soit inferieure ou egale a la somme de toutes 
les valeurs propres diminuee du bruit additif. 

15 De plus, 1 1 estimation du bruit additif est r6alis6e en 

normalisant le bruit instantane qui est evalu6 au moyen du 
signal de reception, du signal de reference et d'une 
estimation de la reponse impulsionnelle du canal de 
transmission. 

20 Avantageusement, en notant A la matrice d 

transmission associee au signal de reference, l 1 expression 
du bruit instantane est la suivante : N 0 = S - A.X. 

Par ailleurs, quelle que soit 1" option eventuellement 
retenue, la matrice de caracterisation r^sulte d f un 
25 operation de lissage. 

Selon un mode de realisation pref erentiel , le signal 
de caract6risation est une estimation de la repons 
impulsionnelle du canal de transmission. 

On peut egalement prevoir que le signal d 
30 caracterisation soit le signal de reception. 

La presente invention apparaitra maintenant de inanier 
plus detainee dans le cadre de la description qui suit od 
sont proposes des exemples de mise en oeuvre £ titre 
illustratif , c ci en r6f6r nee aux figur s ann x6 s qui 
35 representent : 
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- la figur 1, une premiere variante de realisation d 
1 1 invention , et 

- la figure 2, une deuxifeme variante. 

Les elements communs aux deux figures sont affectees 
5 d'une seule et meme reference. 

Le recepteur a d6ja acquis une synchronisation 
grossi&re sur le signal de reception, de I'ordre du temps 
bit f au moyen de I 1 une quelconque des solutions disponibles. 

Ce signal de reception correspond h un signal de 
10 reference produit par l'emetteur et connu du recepteur. Ce 
signal de reference peut etre connu & priori, c'est-a-dire 
qu'il s'agit d'une sequence d'apprentissage form6e d 
symboles identifies. II peut egalement etre connu a 
posteriori au moyen de techniques generiquement reference s 
15 sous le terme de sondage aveugle. Dans ce cas, au cours de 
la procedure de synchronisation, le recepteur regenfere la 
suite des symboles formant le signal de reference & partir 
du signal de reception. 

II convient en premier lieu de caract6riser la 
20 transmission entre l'emetteur et le recepteur. A cet effet, 
on seiectionne un signal source qui, produit par I'emetteur, 
donne au niveau du r6cepteur, aprfes transmission dans 1 
canal, un signal de caracterisation. 

Naturellement , si le signal source est le signal de 
25 reference, ce signal de caracterisation est le signal de 
reception lui-meme. Ce n'est pas cependant pas tou jours la 
solution optimale quant & la complexite et aux performances 
du precede de 1' invention. 

Une autre solution consiste k retenir une impulsion 
30 moduiee comme signal source, le signal de caracterisation 
devenant alors la reponse impulsionnelle du canal de 
transmission . 

A titre d^xemple, le systeme de radiocommunication 
c llulaire numerique GSM fait appel a une s6qu nee 
35 d'apprentissag de 26 symboles, la r6ponse impulsionnelle 
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etant g6n6ralement estim'e avec 5 coefficients puisqu'on 
admet que la dispersion du canal vaut 4. 

Dans ce cas, le signal de reception a une dimension 
maximale de 22 qui est sensiblement plus importante que 
5 celle de la r6ponse impulsionnelle. 

On examinera done successivement deux variantes de 
mise en oeuvre de 1* invention en commengant par le cas oil le 
signal source est une impulsion modulee en reference a la 
figure 1. 

10 L' estimation de la r6ponse impulsionnelle ne pose pas 

de difficult6s en soi car de nombreuses methodes permettent 
d'y parvenir, par exemple la m^thode dite du critfere des 
moindres car r 6s qui est d6crite notamment dans les demandes 
de brevet FR 2696604 et EP 0564849, En matifere de rappel, 

15 cette technique fait appel S une matrice de mesure A 
construite a partir de la sequence d'apprentissage TS de 
longueur n. Cette matrice comprend (n-d) lignes et (d+1) 
colonnes, d representant la dispersion du canal. L 1 element 
figurant h la ifeme ligne et & la j&me colonne est le (d+i- 

20 j)ifeme symbole de la sequence d'apprentissage, soit en 
notant a^ le ieme symbole d'une sequence TS de 26 symboles : 



A = 



f a4 a3 a2 ai ao 
as a4 a3 a2 ai 

*6 a 5 a 4 a 3 a 2 
a 7 



^25 a 21 ) 

La sequence d'apprentissage est choisie telle que la 
25 matrice A fc A soit inversible oil l'op§rateur . fc repr6sente la 
transposition . 

Dans le signal de reception S, on ne prend pas n 
compt les quatr premiers symbol s sq & S3 car ceux-ci 
d6p ndent egalement de symbol inconnus 6mis avant la 
30 s§qu nee d 1 apprentissage , 6tant donn§ que la disp rsion du 
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canal vaut 4. Par un abus de langage on d6finira done 
dorSnavant 1 signal de reception coirane un vect ur S ayant 
pour composantes les symbol es regus, S4, s^, s 6/ ... , s 2 5. 

Dhs lors, 1 • estimation de la reponse impulsionnelle X 
5 prend la forme suivante : 

X « (A* A)* 1 A fc . s 
L' 6 tape suivante du proc6d6 de 1' invention consiste a 
6tablir une statistigue de cette reponse impulsionnelle. Par 
statistique, on entend un ensemble de donnees refl#tant la 
10 valeur moyenne de cette reponse sur une p&riode d 1 analyse. 

On construit done une matrice de lissage des 
diff6rentes estimations X obtenues pendant la period 
d* analyse pour obtenir une estimation de la covariance 
associ£e & cette reponse impulsionnelle. On entend ici 
15 lissage dans un sens tres g£n€ral, e'est-a-dire toute 
operation permettant de lisser ou de moyenner la reponse 
impulsionnelle sur la pSriode d 1 analyse. 

Un premier exemple de lissage consiste 4 effectuer la 
moyenne de la matrice XX* 1 sur la p6riode d 1 analyse suppos6e 
20 comprendre m sequences d 1 apprentissage : 

L(XX h ) = ^ ^ XX h 

L 1 opera teur .h represente la transformation 
hermitienne ou transposition conjugaison complexe. 
25 Un second exemple de lissage consiste k actualiser, 

apres reception de la i&me sequence d 1 apprentissage, la 
matrice de lissage Li«i(XX h ) obtenue a la (i-l)ieme sequence 
d 1 apprentissage au moyen d'un coefficient multiplicatif a, 
ce facteur 6tant generalement connu sous le nom de facteur 
30 d'oubli de lissage et 6tant compris entre 0 et 1 : 

Li(XX h ) = aXiX ± h + (l-a)Li.!(XX h ) 
L f initialisation peut se faire par tous moyens, 
notamment au moyen de la premidre estimation X obtenue ou 
bi n par une moy nne obtenu comm ci-d ssus pour un faible 
35 nombre de sequences d 1 apprentissage. 
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Dans un souci de simplification, la matric d lissage 
L(XX h ) qui est en fait une matrice de caracterisation 
statist ique, sera desormais notee L. 

Le proced6 comprend ensuite une etape de recherche d s 
5 couples (valeur propre, vecteur propre) de la matrice de 
caracterisation . 

Cette etape ne sera pas plus detaillee car bien connue 
de l'homme de metier. 

Les valeurs propres sont maintenant classees , par 
10 ordre d6croissant. En effet, la somroe de ces valeurs 
correspond a 1 ' energie du signal de caracterisation X 
composee pour partie d'un signal utile qui est 1' image du 
signal source et pour partie du bruit additif N du canal de 
transmission . 

15 il vient que les valeurs propres dominantes, celles 

qui sont le plus eievees, representent le signal utile, 
tandis que les valeurs propres les plus faibles representent 
le bruit. 

Selon une premiere option, le proc6d6 consiste h 
20 retenir un nombre predetermine de valeurs propres 
dominantes. Par exemple, pour une reponse impulsionnelle a 
5 coefficients, on retient les deux premieres valeurs 
propres \\ et 

Selon une deuxifeme option, on considfere que le signal 
25 utile pr6sente une 6nergie qui est une fraction 
pr6d6terminee f de l 1 energie du signal de caracterisation. 
Ainsi en notant les valeurs propres pour i variant de 1 a 
p, il y aura d valeurs propres dominantes, d 6tant obtenu 
comme suit : 

30 

d d + 1 

Z fc, Z ^ 

i = 1 ^ ^ . i = 1 ^ 
__ s f et — > f 

Z X. Z K 

i=l 1 i=l 

La fraction f peut §tre fix6e a priori, a une val ur 

de 95 % par exemple. Cett fraction p ut egalem nt deriv r 
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du rapport signal a bruit du signal de reception obt nu par 
ailleurs. 

Selon une troisifeme option du proc€d£, sans doute la 
plus performante, le bruit additif N est estiroS directement 
5 & partir du signal de reception et de la matrice de mesure 
A. En effet, en notant Nq le vecteur bruit affectant le 
signal de reception, il vient que : 
S = AX + N 0 

Compte tenu du fait que les vecteurs S et Nq ont 2 2 
10 composantes, le bruit additif N peut s'exprimer de la 
manifere suivante : 

N = (^) (S - AX)h ( S - AX) 

Naturellement, cette estimation du bruit additif peut 
etre moyennee en liss£e. 
15 En reprenant les notations prdcedentes et n 

normalisant les energies, il vient que : 

P P 
Z *s £ N et Z X. < N 

l = d i = d + 1 

On obtient done les valeurs propres dominantes & 

20 partir d'une estimation directe du bruit. 

Quelle que soit 1' option precedemment retenue, l'6tape 
suivante du proc€d6 consiste & calculer la fonction de 
correlation du signal source avec la somme des vecteurs 
propres associes aux valeurs propres dominantes X^. 

25 Le signal source est sur€chantillone par rapport au 

temps bit et on le notera done g(t) ou t qui repr6sente le 
temps est une variable discrete dont le pas de 
quantification vaut, & titre d'exemple, 1/32 temps bit. II 
est repr€sent£ par un vecteur de meme dimension que le 

30 signal de caract6risation, soit 5 dans l'exemple adopts. La 
fonction de correlation c(t) est calcul6e par exemple pour t 
variant de -1 fit +1 temps bit au moyen de l 1 expression 
suivante : 
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c(t) = z g(t). v. 
i = l 

Le point figurant entre le signal source g(t) et le 
vecteur propre represente classiquement le produit 

scalaire. 

5 La derniere 6tape du proced6 consiste a rechercher la 

valeur tg de t la plus proche de zero qui correspond au 
premier maximum relatif de la fonction correlation c(t) . 
C'est cette valeur particuliere to qui donne l'ecart d 
synchronisation recherchee par rapport au signal de 
10 reception. 

Par ailleurs, on peut considerer la fonction 
complementaire c'(t) suivante : 

c'(t) = 2 g(t). v. 
i = d + 1 

II faut remarquer que la valeur particulifere to 
15 mentionnSe ci-dessus peut egalement §tre obtenue n 
recherchant la valeur de t la plus proche de zero qui 
correspond au premier minimum relatif de la formation 
complementaire c'(t). 

Ces deux methodes pour obtenir l'^cart d 
20 synchronisation to sont done 6quivalentes . 

En r^feence a la figure 2, consid§rons maintenant un 
deuxifeme variante de 1' invention selon laquelle le signal de 
caract6risation est le signal de reception S, si bien que 1 
signal source est maintenant le signal de reference , soit 
25 dans le cas du GSM, la sequence d'apprentissage TS moduli 
GMSK (pour "Gaussian Minimum Shift Keying"). 

La statistique du signal de caract6risation est done 
estimee au moyen d'une matrice de caracterisation qui est 
maintenant obtenue par lissage des differentes occurrences 
30 du signal de reception S. Ici encore, on considere le terme 
lissage dans un sens trfes general. 

La matrice de caracterisation L prend done la forme 
suivante : 
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v. 1 m 
L(SS n ) = — V SS ft ou bien, 
m *f 

Li(SS h ) = aSiSi h + (l-a)Li. 1 (SS h ) 
On recherche ensuite les p f valeurs propres X* i de 
cette matrice et, comrae dans la premiere variante, on 
5 identifie les d 1 valeurs propres dominantes. 

On calcule maintenant la fonction de correlation f(t) 
de Id secjuencs d f apprsntissage iriCd.nl so Gt dc lc s cirjr»c dss 
vecteurs propres v'^ associes aux valeurs propres dominantes 
A* ' ^ • 

10 La encore, la sequence d'apprentissage g 1 (t) est 

surechantillonnee et elle est representee par un vecteur de 
22 composantes. La fonction de correlation devient done : 

d ' 

f(t) = z g f <t).v. 

i = 1 1 

15 Comme auparavant, on peut definir une nouvell 

fonction complementaire f f (t) : 

P' 

f(t) = z g^tj.v. 

i = d' + 1 

Le proc6de se termine de la m§me manifere en 
20 recherchant le premier maximum de la fonction de correlation 
f(t) ou le premier minimum de la fonction complementaire 

t r<t). 

L' invention peut ainsi etre mise en oeuvre de 
differentes maniferes, le point essentiel etant de disposer 
25 d'un signal source et du r6sultat de sa transmission, & 
savoir du signal de caracterisation . 
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REVEND I CAT I ON S 

l)Procede de synchronisation fine sur un signal de reception 
(S) correspondant a un signal de reference (TS) 6mis dans un 
canal de transmission , caracterise en ce qu'il comprend les 
5 etapes suivantes : 

- selection d'un signal source produisant un signal de 
caracterisation (X, S) suite a son passage dans ledit canal 
de transmission, 

- etablissement d'une matrice de caracterisation (L) pour 
10 estimer la covariance dudit signal de caracterisation (X, 

S), 

- identification des valeurs propres dominantes qui sont les 
valeurs propres (X^, X'j.) les plus elev6es de cette matrice 
de caracterisation (L), 

15 - calcul de la fonction de correlation (c(t), f(t)) dudit 
signal source avec la somme des vecteurs propres (vi, v f i) 
associes auxdites valeurs propres dominantes, 

- recherche du premier maximum de cette fonction de 
correlation (c(t), f(t)). 

20 2)Procede selon la revendication 1, caracterise en ce que 1 
nombre (d f d 1 ) desdites valeurs propres dominantes (Xjj, \* ±) 
est predetermine, 

3) Proc6de selon la revendication 1 caracterise en ce que le 
rapport de la somme desdites valeurs propres dominantes a la 

25 somme de toutes les valeurs propres est sup6rieur ou egal & 
un nombre predetermine. 

4) Proc6d6 selon la revendication 1 caracterise en ce que, 
comprenant de plus une etape d' estimation du bruit additif 
(N) dans le canal de transmission, lesdites valeurs propr s 

30 dominantes sont telles que leur somme soit inferieure ou 
egale a la somme de toutes les valeurs propres diminuee 
dudit bruit additif (N) . 

5) Proced6 selon la revendication 4, caracterise en ce que 
l 1 estimation du bruit additif (N) est r6alis6e en 

35 normalisant le bruit instantane (N 0 ) qui est 6valu6 au moyen 
dudit signal de reception (S), dudit signal de reference 
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(TS) et d'une estimation de la reponse impulsionnelle (X) du 
canal de transmission* 

6) Proc6d6 selon la revendication 5 caracteris6 en ce que, en 
notant A la matrice de transmission associee audit signal de 

5 reference (TS), 1* expression du bruit instantan^ (Nq) est la 
suivante : Nq = S - A.X, 

7) Procede selon la revendication 6, caracteris6 en ce que 
ledit bruit additif (N) est de plus moyenn6. 

8) Proc6de selon l'une quelconque des revendications 1 a 7, 
10 caracterise en ce que ladite matrice de caracterisation (L) 

r6sulte d'une operation de lissage. 

9) Proced6 selon l'une quelconque des revendications 
precedentes caracterise en ce que ledit signal de 
caracterisation est une estimation de la reponse 

15 impulsionnelle (X) du canal de transmission, 

10) Proc#de selon l'une quelconque des revendications 1 & 8, 
caracterise en ce que ledit signal de caracterisation est 
ledit signal de reception (S). 
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g(t) 
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c(t) = 2 g(t). v. 
i = 1 



Max[c(t)] -» t 0 
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X«i i = 1, . . . , p 



v'i i = 1 , . . . , d 



f(t) = Z g'(t).V. 
i = l 



Max[f(t)] -> t 0 



